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The two-neutrino double-/? Gamow-Teller and Fermi transitions are studied within an exactly 
solvable model, which allows a violation of both spin-isospin SU(4) and isospin SU(2) symmetries, 
and is expressed with generators of the SO(8) group. It is found that this model reproduces the 
main features of realistic calculation within the quasiparticle random-phase approximation with 
isospin symmetry restoration concerning the dependence of the two-neutrino double-/? decay matrix 
elements on isovector and isoscalar particle-particle interactions. By using perturbation theory an 
explicit dependence of the two-neutrino double-/? decay matrix elements on the like-nucleon pairing, 
particle-particle T — 0 and T — 1, and particle-hole proton-neutron interactions is obtained. It is 
found that double-/? decay matrix elements do not depend on the mean field part of Hamiltonian and 
that they are governed by a weak violation of both SU(2) and SU(4) symmetries by the particle- 
particle interaction of Hamiltonian. It is pointed out that there is a dominance of two-neutrino 
double-/? decay transition through a single state of intermediate nucleus. The energy position of 
this state relative to energies of initial and final ground states is given by a combination of strengths 
of residual interactions. Further, energy-weighted Fermi and Gamow-Teller sum rules connecting 
AZ = 2 nuclei are discussed. It is proposed that these sum rules can be used to study the residual 
interactions of the nuclear Hamiltonian, which are relevant for charge-changing nuclear transitions. 


I. INTRODUCTION 

With increasing sensitivity of double-/? decay {(3P) ex¬ 
periments looking for a signal of Majorana neutrino mass 
the problem of reliable calculation of neutrinoless double¬ 
beta decay (0://?/?-decay) matrix elements becomes 
more urgent [T] . As far as is known their value can not be 
related to any observable and must be calculated by us¬ 
ing tools of nuclear structure theory. Many sophisticated 
nuclear structure approaches including the large basis in¬ 
teracting shell model Hi, the interacting boson model 
d], the projected Hartree-Fock-Bogoliubov method [5], 
the energy density functional method [B] and various ver¬ 
sions of the quasiparticle random phase approximation 
were used to calculate them. The difference among 
obtained results are at the level of factor 2-3 for particular 
nuclear systems [1] . They can be attributed to truncation 
of the nuclear Hamiltonian, many-body approximations, 
and various sizes of the single-particle model space. 

The importance of the improvement of the calculation 
of the 0^'/?/?-decay nuclear matrix elements is accepted 
worldwide. The quality of nuclear structure models can 
be improved by complementary experimental informa¬ 
tion from related processes like two-neutrino double-/? 
decay (2jz/?/? decay), charge- and double-charge-exchange 
reactions, particle transfer reactions, muon capture, etc. 

The 2 j^/?/? decay [T1 [TUIITT] . 

(A, Z) —i (A, 2^ -|- 2) -|- 2e -I- 2r/e, (1) 

is a process fully consistent with the standard model of 
electroweak interaction. So far it has been observed in 
twelve even-even nuclides in which single-/? decay is ener¬ 
getically forbidden or strongly suppressed m- The mea¬ 


surement of 2z//?/?-decay rates gives us information about 
the product of fourth power of axial-vector coupling con¬ 
stant gA and the squared 2i//?/?-decay matrix element 
which is a superposition of double Gamow-Teller 
(GT) and double Fermi (F) matrix elements, 

^ MjS. (2) 

Here, gv is the vector coupling constant. 

The observed values of are used to study the nu¬ 
clear structure and nuclear interactions associated with 
the Oi^/?/? decays. The calculation of requires a con¬ 
struction of wave functions of the even-even initial and 
final nuclei and of a complete set of J+ = O’*", I"*" states 
in intermediate odd-odd nucleus within a nuclear model. 
These wave functions enter also in the evaluation of the 
neutrinoless double-/? decay matrix elements, which has 
a different form. The problem of a reliable calculation 
of is still not solved. Essentially, calculations per¬ 
formed for nuclei of experimental interest overestimate 
the 2://?/?-decay rate HU- The shell model, which de¬ 
scribes qualitatively well energy spectra, does reproduce 
experimental values of only by consideration of sig¬ 
nificant quenching of the GT operator, typically by 60 to 
70% 0. 

In most quasiparticle random phase approximation 
(QRPA) calculations of the particle-particle inter¬ 
action is adjusted so that the 2i//?/?-decay half-life is cor¬ 
rectly reproduced As a result values be¬ 

come essentially independent of the differences in model 
space, nucleon-nucleon interaction, and refinements of 
the QRPA method. Recently, a partial restoration of 
isospin symmetry was achieved within the QRPA Hi 
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by separating the particle-particle neutron-proton inter¬ 
action into its isovector and isoscalar parts, and renor¬ 
malizing them each separately. The isoscalar strength 
parameter is fit as before to 2iy(3/3-decay rates unlike 
the isovector parameter gp/y^, which is determined by the 
requirement that Mp" = 0 dictated by the isospin sym¬ 
metry of the nucleon-nucleon force. As a consequence, es¬ 
sentially no new parameter is introduced as the strength 
of isovector particle-particle force is close to the pairing 
force. 


The Fermi and GT operators are generators of isospin 
SU(2) and spin-isospin SU(4) multiplet symmetries, re¬ 
spectively. In the case of both symmetries being exact 
in nuclei, the 2vj3j3 decay would be forbidden as ground 
states of initial and final nuclei would belong to differ¬ 
ent multiplets. The isospin is known to be a good ap¬ 
proximation in nuclei. Thus, it is assumed that double 
Fermi matrix element is negligibly small and the main 
contribution is given by the double GT matrix element. 
In heavy nuclei the SU(4) symmetry is strongly broken 
by the spin-orbit splitting. But values of M'qj, deduced 
from the observed 2i//I/?-decay rates are especially small 
for nuclei with large A. It is worth noting that the 2v/3f3- 
decay transition to ground state of final nucleus exhausts 
only about 10“"^ of the double GT sum rule [T3]. The 
existence of an (approximate) underlying symmetry re¬ 
sponsible for the suppression of the 2vj3j3 decay, which is 
assumed to be the SU(4) symmetry, justifies approaches 
based on the perturbative breaking of this symmetry for 
construction of wave functions of nuclear states partici¬ 
pating in double-,^ decay transitions. To this category of 
methods belong the phenomenological approach of Ref. 
m and various versions of the proton-neutron QRPA. 


Whether a discussed behavior of M'qj, is a special prop¬ 
erty of nuclei or just an artifact of the QRPA was dis¬ 
cussed within a schematic model which can be solved 
exactly and contains most of the qualitative features of 
a realistic description m- The vanishing of Mqj was 
identified with a dynamical SU(4) symmetry of Hamilto¬ 
nian. Later this model was exploited to examine isovector 
and isoscalar proton-neutron correlations in the case of 
GT strength and double-/? decay m- In this paper we 
extend this schematic model to allow a violation of both 
spin-isospin SU(4) and isospin SU(2) symmetries. The 
main issue is to discuss explicit dependence of both M'qj 
and Mly on the mean field and different components of 
residual interaction by taking advantage of perturbation 
theory. We note that a similar study, which has been 
found to be very instructive, was performed for M'jy by 
discussing violation of isospin symmetry of Hamiltonian 
expressed with generators of the SO(5) group [T7] . 


II. 2://3d-DECAY RATE AND THE 
IMPORTANCE OF THE ENERGY 
DENOMINATORS 


The 2i//?/3-decay occurs as a second-order perturba¬ 
tion of the weak interaction within the minimum stan¬ 
dard model independently of whether neutrinos are Dirac 
or Majorana. The effect of neutrino mixing and masses 
can be safely neglected. The most favorable is the two- 
nucleon mechanism where the successive /? decays of two 
neutrons in the even-even nucleus trigger the 2i//3/? decay. 

The inverse half-life of the 2vf3/3-decay transition to the 
O'*' ground state of the final nucleus is given as follows: 




me 

Stt'^ In 2 


{G.miri^^ (0+) , 


(3) 


where = Gp cos 9c [Gp is Fermi constant and Oq is 
the Cabbibo angle), me is the mass of electron, and 
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f* Ei — Ef — Ee 2^ — E, 


ElElA^''dE,p (4) 


Here, E^,^ = Ei — Ef — Ee^ — Flea ~ E^^ due to energy 
conservation. Ei, Ef, Ee- [Ee, = ^JpI^ + ml) and E„. 
{i = 1,2) are the energies of initial and final nuclei, elec¬ 
trons and antineutrinos, respectively. F(Zf,Ee,) denotes 
relativistic Fermi function and Zf = Z -\-2. A'^’' consists 
of products of nuclear matrix elements, which depend on 
lepton energies: 
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Fn = (0+ II II 0+)(0+ II ^r- II 0+), 

m m 
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and energy denominators are 


Kn{J+) 


Ln{J+) 


2 


{2En{J+)-E,-Ef) + eK 
2 

(2S„(J+) -E,- Ef) - eK 
2 

{2E^{J+)-E,-Ef) + eL 
2 

{2Er,iJ+)-E,-Ef)-eL' 


Here, |0+), |0j^) are the 0+ ground states of the initial 
and final even-even nuclei, respectively, and |0+) (|1^)) 
are all possible states of the intermediate nucleus with 
angular momentum and parity J'^ = O'*" (I’’") and ener¬ 
gies iln(0+) (£;„(!+)). CK = Ee^ + E^^ - Eei - E^^ and 
El = Eei + -Ei /2 — Ee^ — -El /1 ■ We note that formally in 
the limit 2E„ — E^ — Ef =0 one ends up with = 0. 
The maximal value of \eK\ and IclI is the Q value of the 
process. For 2^/3/3 decay with energetically forbidden 
transition to intermediate nucleus {En — Ei > —me) the 
quantity 2E„(J+) — Ei — Ef = Q + 2me + 2(E„ — Ei) is 
always larger than the Q value. We clarify later that this 
quantity can be expressed as a combination of residual 
interactions of nuclear Hamiltonian. 

The calculation of the decay probability is usually sim¬ 
plified by an approximation 


E„( J+) ~ L„( J+) 


2 

Er,{J+)-{E,+Ef)/2- 


( 8 ) 


Then we obtain 


= (9) 

with the Fermi and GT matrix elements given by 

„2i/ ^(o;i|7^-||o+)(o+||r-||o+) 

^ 4^ E„(0+)-(Ei + E;)/2 ’ 

V- (0/ II || ln)(ln II Em|| 0^ 

E„(l+)-(Ei + E/)/2 

( 10 ) 


Here, T~ = 'bn i® the total isospin-lowering opera¬ 
tor. As a result of the above approximation, the separa¬ 
tion of phase space factor and nuclear matrix elements is 
achieved. 

The calculation of M'p' and needs to evaluate 

explicitly the matrix elements to and from the individual 
|0+) and |1+) states in the intermediate odd-odd nucleus, 
respectively. In the shell model and IBM calculation of 
these matrix elements the sum over virtual intermedi¬ 
ate nuclear states is completed by closure after replacing 
En{J^) by some average value E„(J+)): 


M'p' = = 




EjO+)-{E, + Ef)/2’ 




Mqt-cI 


E„{l+)-{E, + Ef)/2 


( 11 ) 


with 


= (o+|r-r-|o+), 

= (0/ I E ' a„|0+). (12) 

m,n 


The validity of the closure approximation is as good as 
the guess about the average energy to be used. This 
approximation might be justified in the case where there 
is a dominance of transition through a single state of the 
intermediate nucleus. 

The T~ operator connects states only in the same 
isospin multiplet. is non-zero only to that extent 

that Coulomb interaction mixes states of different multi- 
plets. As an example 2^/3/3-decay transition "‘^Ca —>■ ^®Ti 
can be considered. The ground state of parent and 
daughter nuclei can be identified with T = 4, Mt = 4 
and T = 2, Mt = 2 states, respectively. A crude esti¬ 
mate of the mixing of the T = 2, Mt = 2 state with 
T = 4 Mt = 2 analog of the ^®Ca ground state due to 
the isotensor piece of Coulomb force implies a negligi¬ 
ble small value Mp'^_^i < 0.02 for this and some other 
2j^/3/3-decay transitions nnj. 

The GT operator T~(Jn connects states only within 
the same spin-isospin multiplet of the SU(4) symmetry, 
which leads to new conserved quantum numbers in addi¬ 
tion to those of spin and isospin. The ground state of the 
initial (A, Z) even-even nucleus belongs to the multiplet 
[n, n, 0] with spin 5 = 0 and isospin T = n = {N — Z)/2 
and it is the only state of this nucleus belonging to that 
multiplet. In the neighbor (A, Z + 1) odd-odd nucleus 
there are two states of the multiplet [n, n, 0], namely the 
isobaric analog state with T = n, S = 0 and the GT 
state with T = n — 1, 5 = 1. In the final (A, Z + 2) 
even-even nucleus, the states belonging to the [n, n, 0] 
multiplet are the double isobaric T = n, 5 = 0 and 
two GT states with T = n — 2, S = 0 and T = n — 2, 
5 = 2. The ground state of the final (A, Z + 2) even- 
even nucleus with T = n — 2,5 = 0 belongs to the 
multiplet [n — 2,n — 2,0]. The SU(4) limit results in 
vanishing matrix elements Mq^ and M^p_^i. The non 
zero double GT matrix element requires a breaking of 
the SU(4) symmetry able to mix the ground and ex¬ 
cited states of the final nucleus. The dynamical origin 
of breaking the SU(4) symmetry is associated with the 
spin-orbit and the tensor potentials which affect mainly 
the mean field. Another possibility is the difference be¬ 
tween strength triplet-singlet, triplet-triplet, and singlet- 
singlet channels of the central potential. We show later 
that the 2i^/I/3-decay NMEs does not depend explicitly 
on the mean field part of the nuclear Hamiltonian. In 
contrast, mainly the differences between triplet-singlet 
and singlet-triplet (spin-isospin) interactions of nuclear 
Hamiltonian contribute to the 2vf5l3 process. This small 
violation of the SU(4) symmetry will be studied in an 
exactly solvable model within the perturbation theory. 
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III. SCHEMATIC HAMILTONIAN EXPRESSED 
WITH GENERATORS OF SO(8) GROUP 

We consider an exactly solvable model m with a set 
of degenerate single-particle orbitals, characterized by I, 
s = \I2, and t = 1/2. The total number of single-particle 
states is n = + !)■ The model is made solvable by 

building a basis entirely from L = 0 operators, i.e., pairs 


of nucleons with spin S' = 0 and isospin T = 1 and with 
S = 1 and T = 0 are allowed. The Hamiltonian of the 
model is an extension of the Hamiltonian introduced in 
Ref. [IS] and possesses the main qualitative features of 
a realistic Hamiltonian relevant to double-/? decay. It 
contains proton and neutron single-particle terms and 
the two-body residual interaction, which components are 
isovector spin-0, isoscalar spin-I pairing and the particle- 
hole force in the T = I, S = 1 channel. We have 


H — T CpAp Qpair I ^ ^ Al,{Ms,0)A,,o{Ms,0)]+9p,Y.<bEa,t 


a.b 


Ho 


{gpair — 9pp°) ^ A\ Q{Ms,0)Ai^Q{Ms,0) + igpair—gpp ^)^0,l(Oi 0)^O.l(0i 0) ■ 

Ms=-l,0,l 


(13) 


H, 


Here, gpair, gp^^,gp^^, and gpt denote the strengths of 
the isovector-like nucleon spin-0 pairing {L = 0, S = 
0,T = 1,Mt ± 1), isovector proton-neutron spin-0 pair¬ 
ing {L = 0,S = 0,T = 1,Mt = 0), isoscalar spin-I 
pairing {L = 0,5' = I,T = 0), and particle-hole force 
(L = 0, 5 = 1,T = 1), respectively. The proton num¬ 
ber operator /Vp, neutron number operator A„, particle- 
particle operators A^g rp{Ms, Mt) , and particle-hole GT 
operators Ea^b are defined in Appendix [X| 

The six particle-particle operators Ag Mt) and 

their Hermitian conjugates together with nine particle- 
hole GT operators i?a,h, total spin 5 and isospin T op¬ 
erators, and total particle number operator (defined for 
convention as Qo = \ + ^n)) represent 28 oper¬ 

ators which generate the group SO (8) [T5|. In case of se¬ 
niority zero, which we will henceforth assume, the SO(8) 
irreducible representation is specified by 7 numbers: (i) 
spatial degeneracy number of levels H 2^ + 1; (h) 

eigenvalue of Qo operator A = il — N/2; (iii) n which cor¬ 
responds to the irreducible SU(4) representation [n, n, 0]; 
(iv) total spin number 5; (v) total spin projection Mg; 
(vi) total isospin number T; and (vii) total isospin projec¬ 
tion Mt- As we are constrained by the set of degenerate 
I shells with total degeneracy H and given particle num¬ 
ber N, for basis state we introduce the abbreviation as 
follows: 


|5,Ms,T,MT,n), or |5Tn) . 


(14) 


We note that matrix elements of generators SO(8) group 
are known in this basis [TS] [TUj , which allows diagonaliza- 
tion of the Hamiltonian (131. Relevant expressions can 
be found in Apppendix|B 

The physics associated with a simplified version of the 
Hamiltonian in Eq. (13) was studied previously with 
emphasis on energy levels [HIMIIII], the extreme sensi¬ 


tivity of Mqt to the strength of proton-neutron particle- 
particle interaction m, and the interplay between the 
isoscalar and isovector pairing models [TB]. Here, we dis¬ 
cuss the role of the violation of the isospin SU(2) and 
spin-isospin SU(4) symmetries and of different compo¬ 
nents of Hamiltonian in the calculation of two-neutrino 
double-,0 decay matrix elements by taking the advantage 
of the perturbation theory. 


The Hamiltonian in (13) is decomposed in two parts: 
Hq, the unperturbed Hamiltonian, and 77/, the perturb¬ 
ing one. The eigenstates of unpertubated Hamiltonian 
Hq are characterized by the number of nucleon pairs 
(only systems with even number of nucleons are consid¬ 
ered), the isospin T, and a quantum number n corre¬ 
sponding to the irreducible SU(4) representation [n, n, 0]. 
The possible values of quantum number n for system with 
N particles in the set of degenerate I shells with degene¬ 
racy fl and given T, 5 are 5 -I- T, 5 -I- T -|- 2, • • •, nj^ax^ 
where Umax = M/2 if A/’/2 < H and Umax = 2H — M12 
otherwise [16]. The single-particle and particle-hole in¬ 
teraction components of 77o violate both isospin and 
spin-isospin symmetries and as a consequence energies 
of states with the same quantum numbers T and n are 
different for a given = {N — Z)/2 (T^ = Mt)- N and 
Z are numbers of neutrons and protons, respectively. If 
gpair = and gpair = gpp^ the isospin and spin- 

isospin symmetries of particle-particle interaction are re¬ 
stored we get H = Hq and 77/ = 0. If gpair ^ gp^° 
and/or gpair ^ gpE\ the Hamiltonian (13) is not more 
diagonal in basis (141 and states with ditterent quantum 
numbers T and n are mixed. The eigenstate of the Hamil¬ 
tonian \S,Ms^T ’can be expressed with eigen¬ 
states of unperturbated Hamiltonian 77o as follows: 


|5,Ms,r,MT,n') = 
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Sn,T ^,Ms,T,MT,n ^ ^S, T, Mt, n) . (15) 

Here, we assume a small violation of the SU(4) symme¬ 
try, which can be treated by a perturbation theory. The 
prime symbol by quantum numbers T and n {T' and n') 
indicates that these quantum numbers are not more good 
quantum numbers due to the violation of SU (4) symme¬ 
try and that the dominant component in the expansion 
over states with a good isospin and the SU(4) quantum 
number is that with T' = T and n' = n. A diagonal- 
ization of Hamiltonian requires calculation of matrix el¬ 
ements 

(5, Ms, T, Mt, n±2\H\S, Mg, T, Mt, n ), 

(5, Ms, T±2, Mt, n\ H\S, Mg, T, Mt, n), 

{S, Mg, T±2,MT,n±2\H\S, Mg, T, Mt, n). 

The corresponding expressions are given explicitly in 
Apppendix 

We shall assume a small violation of the SU(4) sym¬ 
metry due to Hi ^ 0, namely Qpair — and/or 

gpair — 9p^^ ■ For the numerical application we consider 
a set of parameters as follows m- 

ep = 1.2MeV e„ = l.lMeV H = 12, 

J\f — 20, 9pair — O.SMeV, gph — 1.5 gpair' 

(16) 

The initial, intermediate and final states of the double¬ 
beta decay transition will be identified with the isospin 
projection = 4, 3, and 2. For these three values of 
Tj the corresponding numbers of neutrons and protons 
{N,Z) are (14,6), (15,7), and (12,8), respectively. 

In Fig. [^we present states with energy EgTn of dif¬ 
ferent isotopes. The considered level scheme illustrates 
the situation with double GT transition for ^®Ca as the 
isospin and its projection of the initial and final ground 
states correspond to those of ^®Ca and "^®Ti. We note that 
in nuclear physics the isospin symmetry is conserved to 
a great extent. Within the studied model in the SU(4) 
symmetry limit the ground states of “^^Ca and ^®Ti can 
be identified with S = 0, T = 4, = i, n = 4, and 

S = 0, T — 2, Tz = 2, n = 2, respectively, and the inter¬ 
mediate states in ^®Sc with S=1 (T=3,5,7, and 9) T2=3 
(n=4, 6, 8 and 10). As the GT operator is not chang¬ 
ing quantum number n, the double GT matrix elements 
connecting initial and final ground states is nonzero only 
to the extent the breaking of SU (4) symmetry mixes the 
high-lying (0,4,4) analog of the "^^Ga ground state into 
(0,2,2) analog of the ^STi. 


IV. DOUBLE FERMI AND GT MATRIX 
ELEMENTS WITHIN THE PERTURBATION 
THEORY 


We study the double GT and Fermi matrix elements 
using perturbation theory within the discussed model 
close to a point of restoration of the SU(4) symmetry of 



T 

Z 


FIG. 1. (Color online) Eigenenergies EsTn of the Hamilto¬ 
nian (131 for set of parameters (|16|l, Tz=4, 3, and 2 and by 


assuming Hi = 0. Energy states are labeled by spin, isospin, 
and the SU(4) quantum number n: {S,T,n). The levels with 
different value of S' + T are displayed in different colors on¬ 
line: S-l-r = 2 (green), 4 (red), 6 (orange), 8 (black), and 10 
(blue). 


particle-particle interaction of H. First, we assume a con¬ 
servation of the isospin symmetry by the particle-particle 
interaction and a subject of interest will be M'qj. as func¬ 
tion of the isospin of the initial state. Then a weak vio¬ 
lation of the isospin symmetry is allowed and the depen¬ 
dence of M]Z and M'qj, on both quantities gpair — 
and gpair — 9p^^, which violates the SU(4) symmetry, is 
analyzed. 


A. The GT matrix element in the case of isospin 
symmetry 

We consider a small violation of the SU(4) spin-isospin 
symmetry in nuclear Hamiltonian ( |13[ ) due to gpair 7 ^ 
gp^^ and that isospin is a good quantum number, i.e., 
9pazr = which implies Mjy = 0. 

As an example we discuss in details the GT matrix 
element for 2i'j3f3 decay from the state with S = 0, T = 
Mt = 4 to the state with S = 0, T = Mt = 2. The 
corresponding transition is 

|0, 0,4,4,4) ^ |1, Mg, 3, 3,4) ^ |0,0, 2, 2, 2) . (17) 

In the case gpair = 9p^^ one finds that M^l/ = 0 as 
eigenstates of the GT operators are diagonal in SU(4) 
quantum number n and the initial and final states are 
assigned into different SU(4) multiplets. By breaking the 
SU(4) symmetry of particle-particle interaction the quan¬ 
tum number n is not more a good quantum number and 
states with different n are mixed. By keeping in mind a 
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small violation of this symmetry we denote perturbated 
states and their energies with a superscript prime sym¬ 
bol {\S,Ms,T,MT,n'), i^s,Ms,T,MT.n)> unlike the states 
with a definite quantum number n {\S, Ms,T, MT,n), 
Es,Ms,T,MT,n)- 


Up to the second order of parameter {gpair — 9p^^) we 
get (for sake of simplicity a shorter notation of states and 
energies (14) is used) 


E'o22 = Eo22 + (0221 Hi |022) + 


(0241 iJ/ 1022) 


i?022 — 1?024 
= -\- 8Cp 94:gpair ~h Qgph 

132 {gpa^r - glp^ f 8316 


+ {gpair gpp ) 


lO^pair 4” ‘^9gph 25 


(18) 


= £^134 + (134| 77, |134) + 


(136|77,|134) 


Ei34 — Ei3q 
— 13e7i Tcp ^Agpair 4“ ^8gph 

201 {gpa^r - 9l=°f 10125 


{gpair gpp ) j 


145p 


28gph 49 


(19) 


E'o44 = Eo44 + (044| Hi |044) 


(046| Hi |044) 


£-044 — £'046 

= 14en 4“ Ocp 8^gpair 4- 12^p/i 


4- {gpair gpp )- 


108 {gpazr - gtp^ f 7425 
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lAgp 


‘28gph 49 


( 20 ) 


For the double GT matrix element we have 



FIG. 2. (Color online) Matrix element Mgy for the double- 
GT two-neutrino double-d decay m ode as function of ratio 
9 pir^/gpair for a set of parameters (16l. Exact results are 
indicated with a solid line. The results obtained within the 
perturbation theory up to the Erst and second order in Hi 
contribution to Hamiltonian (131 are shown with dashed- 
dotted and dashed lines, respectively. The restoration of spin- 
isospin symmetry of particle-particle interaction is achieved 

fcr Qpp jgpair ~ 1. 


denominator in (21) we obtain 


(^£'134 Eq 22 — bg._ 


\ipair 


■ 9gp}i 


{gpair gpp ) - 


39 (5pair-5pp-°)V 1249263 \ 


gpair 4“ 2gpfi y 171500 J 


(23) 




= E 


(022'| ar- |13n') • (13n'| ar" |044') 


£^(3n 


(022'I 


(£'044 £' 022 ) /2 

134') • (134'|CTr- |044') 


£^134 (£^044 £' 022 ) /2 


( 21 ) 


The allowed intermediate states |13n') are those with 
S' = 1, r = 3 and n' = 4,6,8 and 10. We note that 
up to second order of perturbation theory there is only a 
single contribution through the intermediate state |134') 
and the product of two corresponding /3-amplitudes (nu¬ 
merator of (21) takes the form 


( 022 ' 


(TT 


134') • (134' |Sr-|044') = 144W- 

' ' V 35 


{gpair gpp *^) ‘2Q7{gpair gpp 


lO/Zpair 4- 20//p/i 8b{l-0gpair 4- 20//p/^)^ 


( 22 ) 


We see that if gpair =5^ ° GT matrix element vanishes. 
With help of Eqs. (18i, (19) and (20) for the energy 


It is worth noting that neither the numerator nor deno¬ 
minator of depend explicitly on the single-particle 
energies e„ and Cp. If we restrict our consideration to the 
first-order perturbation theory we find 




^'^^\/W(9pair - gpp °) 


{^9p 


9gph){79gpair 4- 2{)gpii) 


(24) 


In Fig. -^GT plotted as function of ratio 
9pp^ /9pair- We see that results obtained within the 
second-order perturbation theory agree well with ex¬ 
act results within a large range of this parameter. For 
9pp^^!9pair = 1 the restoration of the SU(4) symmetry 
of particle-particle interaction is achieved, i.e., Mpyp is 
equal to zero. We notice that if the quantity gp//^^/gpair 
is within the range (0.8,1.2) the first-order perturbation 
theory seems to be sufficient. 

Usually, ground states of stable even-even nuclei are 
identified with isospin T = T^. The dependence of 
on the isospin of the initial nucleus is presented in Fig. 
We see that for a fixed value of gpp'^/gpair (be., break¬ 
ing of the SU(4) symmetry) the absolute value of M'q^ 
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FIG. 3. (Color online) Matrix element for the double- 
GT two-neutrino double-/! decay mode as function of ratio 
/gpair for different initial state with T = Mt {Mt = 2, 
4, 6, 8, and 10). 


FIG. 4. (Color online) Matrix elements Mp" and Mgy as 
function of ratios gpjT^/gpair and g^^^ /gpair for transition 
from the initial |04'4') to final |02'2') ground state and a set 
of parameters ( |16[ ). The results are obtained within the per¬ 
turbation theory up to the second order. 


decreases with increasing isospin T. We note that apart 
from the shell effects of magic nuclei this tendency is 
observed also for measured 2vj3jS-decsy matrix elements 

m- 


B. The Fermi and GT matrix elements in the case 
of broken SU(2) and SU(4) symmetries 


The main task to be addressed in this subsection 
is determining the dependence of and on 

both quantities /Opair and gp^^/gpair- Recall that 

3pp^ 7^ 9pair breaks both the SU(2) isospin and the 
SU(4) spin-isospin symmetries of particle-particle inter¬ 
action unlike gp^^ ^ gpair, which is associated only with 
the violation of the SU(4) symmetry. 

We consider the 2://3/3-decay transition from the initial 
|04'4') to final |02'2') ground state. Up to the first or¬ 
der in the perturbation theory for double Fermi and GT 
matrix elements we find 


Mjy = 

48 ( 

y 5 \ 

gpair gpp ^) 

i^9pair ^9ph){^^9pair ^9ph) 

II 

i44yf 

f {gpair gpp *^) 

^9pair ^9ph 

{\^9pa,ir 


(25) 


-b 


‘^9ph {9pair 9pp ) 


(IbUpair “b “b Qgph) 


(26) 


We see that Mp' depends only on strength of the isovec¬ 
tor interaction gp^^ unlike which depends also on 


the strength of the isoscalar interaction gp^^ ■ Due to the 
violation of the isospin symmetry the final ground state 
|02'2') is mixed with both first |04'4') and second |02'4') 
excited states (see Fig. [^, resulting in gp^^ contribution 

to M^rp. 

In Fig. H we present behavior of M^' and as 

function oigp^^ / gpair {gp^° / gpair) for a particular val¬ 
ues of gpp^/gpair {gpp^/gpair)- Results were obtained 
within the perturbation theory up to the second order. 
We see clearly that for gp^^/gpair = 0 matrix element 
M'P' does not depend on gp^° and varies strongly with 
change of gp^^ ■ A different behavior offers which 

weakly depends on the gp^^ and significantly on the 
gpp^- These conclusions agree qualitatively well with 
results obtained for two-neutrino double-/3 decay transi¬ 
tions within the proton-neutron QRPA with restoration 
of the isospin symmetry [8] . The advantage of the study 
which considered the schematic model and in perturba¬ 
tion theory is that explicit dependence of and M'qj, 
on isoscalar and isovector strength of particle-particle in¬ 
teractions can be determined. 


V. ENERGY-WEIGHTED SUM RULE OF AZ = 2 
NUCLEI 


In Ref. [17] the double Fermi and GT sum rules asso¬ 
ciated with AZ = 2 nuclei were introduced. We have 


oew 

Op 


(*,/) = V(^;„- 




Ef 


){f\T-\n) {n\T-\t) 


2 
























= ^{f\[T-,[H,T-]]\t) (27) 

/) ^ ^(£;„ - (/I dcT In) • (n| Ogt \^) 

n 

= [(C’gt)-m, [H, {Ogt)m]] \i) 

M 

(28) 

with 

A 

Ogt = y^Tfc~'^k, (29) 

fe=i 


where |z > (|/ >) are 0+ ground states of the initial 
(final) even-even nuclei with energy Ei {Ef), and |1+ > 
(|0+ >) are the 1+ (0+) states in the intermediate odd- 
odd nucleus with energies En. 

If there is a dominance of contribution through a sin¬ 
gle or few states of the intermediate nucleus, energy- 
weighted sum rules (27) might be exploited to adjust 
the strengths of the residual interaction of Hamiltonian 
for nuclear structure calculations. The left-hand side of 
Eq. (271 might be determined phenomenologically, unlike 
the right-hand side of Eq. (27), which requires evalua¬ 
tion of the double commutator within a nuclear model 
and can be expressed in terms of the strengths of resid¬ 
ual interaction. Due to a double commutator of nuclear 
Hamiltonian with charge-changing Fermi and GT opera¬ 
tors connecting states with AZ = 2 the energy-weighted 
sum rule S^Qrp(i^ f) does not depend explicitly on the 
mean field part of the nuclear Hamiltonian. 


TABLE I. The coefficients a, b, c and d of t he e xpansion of 

the energy denominator - ^ [see Eq. (341] associated 

with the dominant double GT (double Fermi) transition from 
the initial ground state jOT = Mtu) to the final ground state 
jo T = Mt — 2 n —2) through a single state of the intermediate 
nucleus |1Mt n) (jO Mt n)). Coefficients are presented for 
different isospin T = Mt of the initial state. 


T = Mt 

Transition 


Coefficients 

a 

be d 

2 

GT 

3 

5 -59/15 44/5 


Fermi 

3 

3 50/3 -59/5 

4 

GT 

5 

9 -64/35 39/5 

4 

Fermi 

5 

3 401/35 -192/35 

6 

GT 

7 

13 -71/63 340/63 


Fermi 

7 

3 482/63 -71/21 

8 

GT 

9 

17 -80/99 103/33 


Fermi 

9 

3 469/99 -80/33 

10 

GT 

11 

21 -7/11 12/11 


Fermi 

11 

3 26/11 -21/11 


We analyze Sj^qt for the Hamiltonian (13) and by 
exploiting the commutation relations of the SO(8) group. 
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For the case |i) = |044'), |/) = |022'), we get 
5g"^(04'4',02'2') ^ ^(E(3„ - 1^0444^) 

n 

X (02'2'| ar- |13'n') • (13'n'| ar" |04'4') 

= 6(gJ=° - gpa^r) (02'2'| Ho_i(0, -l)7lo,i(0,1) |04'4') 

-gph (02'2'| ar- ■ ar- |04'4') - Sgph (02'2'| T-T" |04'4') 

(30) 

and 

,5|,-(04'4',02'2') ^ 

71 

X (02'2'|r- |04V) • (04V|T- |04'4') 

= 2[gpa^r - gl=^) (02'2'| -l)4lo,i(0,1) |04'4'). 

(31) 


We note that the dominant contribution to 
S'gf,(044',022') and S'!!"(044', 022') comes from the 
transition through the single intermediate states |134') 
and 1044'), respectively. By exploiting the first-order 
perturbation theory to evaluation of matrix elements 


in Eqs. (31) and (30) for a combination of energies of 


involved states we find 
Eq 44 -I- i?n 




022 


= 


pair 


^gph 


(32) 


64 T=i\ 39 T=o\ 

~ 35 ~ ^PP ^ [gpair ~ gpp ) 


K 


^044 


•Jpp 
^022 


0.0,4,3,4 


= 55; 


pair 


^gph 


(33) 


gg [gpair gpp j gg [gpair gpp )■ 


The result in Eq. (33) is in agreement with above cal¬ 


culated expression for energy denominator in Eq. (23), 


which was derived by assumption of the isospin conser¬ 
vation. 

We see that considered energy-weighted sum rules 
imply useful relations between three energies of states 
appearing in the denominator of the double GT or 
Fermi matrix elements and nucleon-nucleon interactions. 
The energy denominator associated with the dominant 
double-GT (double-Fermi) transition from the initial 
ground state |0T = Mtu) to the final ground state 
|0 T = Mt — 2 n — 2) through a single state of the in¬ 
termediate nucleus |lMr n) (|0 Mt n)) can be written 
as 



(34) 


— P^gpair 4 ^gph 4” C i^gpair gpp ) 4“ d i^gpair gpp ) • 

Here, a, 6, c and d are coefficients. The perturbation 
theory up to the first order is considered. For different 
isospin T = Mt of the initial ground-state coefficients a, 
b, c and d are presented in Table |Tj We see that for larger 
value of T the value of the energy denominator is affected 
less by the violation of both the isospin and spin-isospin 
symmetries as it is for smaller value of T. 
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VI. CONCLUSIONS 

The anatomy of the two-neutrino double-/? decay ma¬ 
trix elements was studied within a schematic model, 
which can be solved exactly and yet contains most of 
the qualitative features of a more realistic description, 
and by taking advantage of the perturbation theory. We 
paid attention to violation of both spin-isospin SU(4) and 
isospin SU(2) symmetries of particle-particle interaction 
of Hamiltonian. The isospin violation originates from the 
difference of the proton-proton and the neutron-neutron 
pairing force compared to the proton-neutron isovector 
pairing force. The break down of the SU(4) symmetry is 
a consequence of the difference of the like-nucleon pair¬ 
ing interaction compared to the proton-neutron isoscalar 
interaction and/or to the proton-neutron isovector inter¬ 
action, which violates also the isospin symmetry. 

By using perturbation theory, an explicit dependence 
of the two-neutrino double-/? decay matrix elements on 
different constituents of the Hamiltonian was established. 
It was found that the mean-held part of the Hamiltonian 
does not enter explicitly in the decomposition of M'p' and 
Mq^p and is related only to the calculation of unpertur- 
bated states of the Hamiltonian. This general conclusion 
is valid for any mean held approximation. In the case 
of medium and heavy heavy nuclei the SU(4) symmetry 
is strongly broken by the spin-orbit splitting, affecting 
strongly the mean held part, unlike the interaction part 
of nuclear Hamiltonian. This fact might be an explana¬ 
tion for a smallness of Mq^ being governed by a small 
violation of the SU(4) symmetry by the particle-particle 
interaction of the Hamiltonian. 

The obtained expressions for Mjy and Mpi/ supported 
by numerical calculation up to the second order in per¬ 


turbation theory conhrm the hnding achieved within 
the proton-neutron QRPA approach that Mp' depends 
strongly on the isovector part of the particle-particle 
neutron-proton interaction, unlike which depends 

strongly on its isoscalar part. By assuming a hxed vi¬ 
olation of the SU(4) symmetry by the particle-particle 
interaction it is shown that value of decreases by 

an increase of the isospin of the initial ground state. This 
tendency is found also in the case of measured double- 
GT matrix elements being partially spoiled by different 
pairing properties. We also showed that contains 

a small component due to violation of the isospin sym¬ 
metry. By keeping in mind that in nuclear physics the 
isospin symmetry is conserved to a great extent it is rec¬ 
ommended for evaluation of double-/? decay matrix ele¬ 
ments to use many-body approaches with a conservation 
or restoration of the isospin symmetry [aiig. 

An important result coming from the analysis within 
perturbation theory is that there is a dominance of 
double-/? decay transition through a single intermediate 
state. Further, the importance of the energy-weighted 
sum rule associated with AZ = 2 nuclei for htting differ¬ 
ent components of the residual interaction of the Hamil¬ 
tonian was pointed out. It goes without saying that fur¬ 
ther studies, in particular those in which a realistic nu¬ 
clear Hamiltonian is used, are of great interest. 
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Appendix A: Operators in SO(8) schematic model 


We consider a set of single-particle states with the associated creation and annihilation operators, 

Oimmsmt) which are labeled by orbital angular momentum I, its projection on z axis m, and z components of spin 
(tos = 1/2) and isospin {rrit = 1/2). 

The particle-particle operators entering the Hamiltonian (13) are given by [T^ 

Ay{M.,.MT)= Y. (S.r) = (0,l)or(l,0), (Al) 

the particle-hole GT operators take the form 


Ea,b = \msmt) 

l,m,ms ,mt 


mrrismt ’ 


and particle number operators are written as 


H ^ TOt„,t^=±l/2. 

Lm,ms ,m,l 


(A2) 


(A3) 
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Here, tXa and r;, represent spherical components of the single-particle Pauli spin and isospin operators with convention 
used in Ref. [18] . 


Appendix B: Relevant SO(8) matrix elements 


The matrix elements of SO(8) operators in the basis of zero-seniority states are given in Refs. [TS] and [121 . Here, 
we present the SO(8) matrix elements relevant for calculation of the double GT and Fermi matrix elements. We have 


(5, Ms, T, Mt, n\ 4.o(A^s, 0)A,,o{Ms, 0) |5, Ms, T, Mr, n) 

(n-|-n-|-A-|-6)(n-n-A) f (S+l)(n-|-S-|-T-|-4)(n+S-T-|-3)-l-S(n-S-|-T-|-3)(n-S-T-|-2) 
8(n-|-2)(n+3) [ (2S-I-1) 


. (O—n-t-A+2) (H-t-n—A-t-4) 

8(F+T)(F+2) 


(S-l-l)(n-S+T-|-l)(n-S-T)+S(n+S-T-|-l)(n+S-l-T+2) 

(2S-I-1) 


{S, Ms, T, Mt, n\ Em. Mt)A^M Mt) \S, Ms, T, Mt, n) 

(n-|-n-|-A-|-6)(n-n-A) 1 (T+l)(n-|-r+S-l-4)(n+T-S-|-3)-l-T(r!,-T-|-S+3)(n-T-S-|-2) 
8(n-|-2)(n-|-3) [ (2T+1) 


. (n-n-|-A-|-2)(n-|-n-A+4) 
+ 8(n+l)(n-|-2) 


(T-|-l)(n-r-|-S-|-l)(n-r-S)+T(n+T-S-|-l)(ra+T-|-S+2) 

(2T-I-1) 


(5, Ms, T, Mt, n + 2\ Em. 4.i(0. ^t)^o.i(0, Mt)\S, Ms, T, Mt, n) 

= - {S, Ms, T, Mt, n + 2\ Em^ A.oi^s, 0)Ai.o(Ms, 0) |^, Ms, T, Mt, n) 

_ 1 / (n+n+A+6) (n+n—A+6)(f2—n+A)(0—n—A) 

8(n+3) Y (n+2)(n+4) 

X + S + T + 4)(n + S-T + 3)(n -S + T + 3)(n -S-T + 2) 


{S, Ms, T, Mt, n\ Ea,b \S, Ms, T, Mt, n) 

= n(n -b 4) - S{S + 1) - T{T + 1) 


(5, Ms, T, Mt, n\ 4 ^(0,0)Ho.i(0,0) |5, Ms, T, Mt, n) 

71 (fi—77-t-A-t-2)(S2-t-7i—A ^4) 

^ 4(F+2) 

'T(n+S+T+2)(n-S+T+l) (T-Mt)(T+Mt) , (T-|-l)(n+S-T-|-l)(n-S-T) (T+M. + l)(r-M.-|-l)' 
2n(n-|-l)(2T-|-l) r(2T-l) 2n(r7-|-l)(2T-|-l) (T+l)(2T-|-3) 

I (n+4)(Q-t-n-t-A-t-6)(Sl —n—A) 

4(77-1-2) 

T(77-S-T+2)(77-|-S-T+3) (T-Mt)(T+Mt) . (T-|-l)(77-S-|-T-|-3)(r7-|-S-|-T+4) (T+M. + 1)(T-M.-|-1) 

2(77-|-3)(77-|-4)(2T-|-1) T(2T-1) + 2(77+3)(77-|-4)(2T-|-l) (T+l)(2T-|-3) 


{S, Ms, T + 2, Mt, n\ AJ ^(0,0)Ao,i(0,0) |5, Ms, TMt, n,) 


_ ! {T+Mt + 1)(T-Mt + 1) 

(2T-|-3)(T-|-1) 


{T+Mt+2)(T-Mt+2) 

(2T+3){T+2) 


■ 77(0-77-1-AH-2)(n-|-77-A+4) / (T-|-2)(77 + S-|-T+4)(77-S-|-T + 3) / (T-|-1)(77 + S-T-|-1)(77-S-T) 

4(77+2) Y 277 ( 77 + 1 ) (2T+5) V 277(77+1)(2T+1) 

(77+4)(0+77+A+6)(0-77-A) /(T + 2)(77-S-T)(77 + S-T+1) / (T +1) (77-S+T + 3) (77+S+T+4) ' 

4(77+2) Y 2(77+3)(77+4)(2r+5) Y 2(77+3)(77+4)(2T+1) 


{S, Ms, T, Mt, n + 2\ 0)Ao,i(0, 0) |^, Mg, T, Mt, n) 


_ j (n+2)(n —n+A)(r2+n—A+6) 

4(n+4) 


(n+4)(r2+n+A+6)(r2—n—A) 
4{n+2) 


X 


T(n+5+T+4)(n-5+T+3) T(n-5-T+2)(n+5-T+3) {T-Mt){T+Mt) 




2(77+2)(77+3)(2T+1) 2(77+3)(77+4)(2T+1) T(2r-1) 

(T+1)(77+S-T+3)(77-5-T+2) (T+1)(77-S+T+3)(77+S+T+4) (T+M. + 1)(T-M. + 1) 


2(77+2)(77+3)(2T+1) 


2(77+3)(77+4)(2r+l) 


(T+l)(2T+3) 
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{S, Ms, T + 2,MT,n + 2| 4 ^(O, 0)Ao,i(0,0) |^, Mg, T, Mr, n) 


_ j (n+2)(n—n+A)(0+n—A+6) 

4(n+4) 


(n4-4) (n+n+A4-6) (S4—n—A) 
4(n+2) 


(T+2)(n+S+T+6)(n-S+T+5) 

2(n+2)(n+3)(2T+5) 


(T+l)(n-S+T+3)(n+S+T+4) 

2(n+3)(n+4)(2T+l) 


(T (T — 447 ’ + 1 ) 

(2T+3)(T+1) 


(T+Mt+2)(T-Mt+2) 

(2T+3)(T+2) 


(0,0, T - 2, Mt - 2, n - 2| 4 ^ (0, -l)Ao.i (0,1) |0, 0, T, Mt, n) 


(w+2)(0+n+A+4)(n-w-A+2) 

An 


—7i+A+2)(0+7i—A+ 4) 
4(n+2) 


(T-l)(n-l+T)(n+T) 

2(n+l)(n+2)(2T-3) 


T(n+T+2)(n+T+l) 

2n(n+l)(2T+l) 


(T+MT-2)(T+Mr-3) 

(2T-2)(2T-1) 


(T+Mr-1)(T+MT) 

(2T-1)2T 


(0,0, T - 2, Mt - 2, n| 4 ^(0, -l)Ao,i(0,1) |0,0, T, Mt, n) 


(T+Mt-2)(T+Mt-3) 

(2T-2)(2T-1) 


(T+Mt-1)(T+Mt) 

(2T-1)2T 


—?i+A+2)(S2+?i—A+4) 
4(n+2) 


T(«+T+2)(«+T+l) 

2n(ra+l)(2T+l) 


(T-l)(n-T+3)(n-T+2) 

2n(n+l)(2T-3) 


I ( 71 + 4 )(£2+71+A+6)(£2— 71 —A) j (T —l)(77+2^+l)(77+7^+2) jT(^n —T^+2)(77—7^+3) 
4(77+2) V 2(77+3)(77+4)(2T-3) V 2 ( 77 + 3 )( 77 + 4 )(2T+1) 


(0,0, T, Mt - 2, n| 4 ^(0, -l)74o.i(0,1) |0,0, T, Mt, n) 


77 (£2—77+A+2)(£2+77—A+4) 
^ 4(77+2) 


T(77+T+2)(77+T+1) /(T-Mt + 1)(T-Mt+2) / (T+Mt-1)(T+Mt) 

2r7(77+l)(2T+l) y 2T(2T-1) Y 2T(2T-1) 


, (T+1)(77-T+1)(77-T) / (T+Mt)(T+Mt-1) / (T-Mt + 1)(T-Mt+2) ' 

2r7(77+l)(2T+l) y (2T+2)(2T+3) y (2T+2)(2T+3) 


(77+4)(£2+77+A+6)(£2—77—A) 
4(77+2) 


T(77-T+2)(77-T+3) /(T-Mt + 1)(T-Mt+2) j (T+Mt-1){T+Mt) 
2(77+3)(77+4)(2T+1) y 2T(2T-1) V 2T(2T-1) 


+ 


2(77+3)(77+4)(2T+1) 

For Tf = T — 2 the GT matrix element is 


(T+1)(77+T+3)(77+T+4) /(T+Mt)(T+Mt-1) / (T-Mt + 1)(T-Mt+2) 


(2T+2)(2T+3) 


(2T+2)(2T+3) 


{0,0,T -2,T -2,n'\aT- \1, Ms,T - 1,T - l,n') ■ {1, Ms,T - 1,T - l,h'\ aM \0,0,T,T,n') 

h' ,Ms 

_ \ '' r) n' h' fi' n' 

= 2_^ ^^0,0,T-2,T-2,n^l,0,T-l,T-l,n^l,0,T-l,T-l,n^0,0,T,T,n 
n,n 

I (T - l)r(h -T + 2)(n -T + 3)(n + T+ l)(h + T + 2) 

(2T-l)(2r+l) 

Within the perturbation theory the subject of interest is the GT matrix element as follows: 

^ (0, 0,4,2,4| (7T- |1, Ms, 3,3,4) • (1, Mg, 3, 3,4| ar" |0, 0,4,4,4) = 

Ms 
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